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Abstract — Consideration is given to attitude control associ- 
ated with capturing a free-flying asteroid using an axisymmetrie 
spacecraft. Asymptotically stable controllers are designed to 
align the spacecraft axis of symmetry with a line of descent 
that is fixed in the asteroid, and to eliminate all relative angular 
velocity before capture takes place. An analytical expression is 
presented for the torque required to maintain alignment of 
the axes of symmetry of the spacecraft and an axisymmetrie 
asteroid. After the asteroid is securely captured, the angular 
velocity of the rigid composite body relative to an inertial frame 
is arrested; we present a controller that is asymptotically stable 
and stays within specified thrust limits. 

I. Introduction 

Capturing and returning an asteroid to the Earth-Moon 
system has garnered a great deal of attention recently (Ref. 
[1]). References [2] and [3] discuss various phases of the 
operations when the spacecraft is in the vicinity of the 
asteroid, including initial approach, hover, spacecraft spin- 
up, descent, ascent (if a small boulder is picked up from a 
large asteroid), and de-spin. Dynamics and control in each 
of these activities is analyzed in order to determine the 
velocity increments and control force and torque that must 
be provided by a Reaction Control System (RCS), and the 
mass of the propellant that will be consumed. 

This paper addresses in depth the attitude control required 
for spin-up and de-spin. First, the spacecraft is reoriented 
and spun up so that it points to the asteroid properly, and 
all angular velocity relative to the asteroid is eliminated to 
ensure capture of the asteroid. Second, the body consisting 
of the spacecraft and the captured asteroid is de-spun to 
eliminate angular velocity relative to an inertial reference 
frame, to prepare for the return to the Earth-Moon system. 

We make the following assumptions. The asteroid is a 
rigid body with mass distributed in the most general possible 
way; that is, it can possess three distinct principal moments 
of inertia, and is therefore referred to as unsymmetric. 
Before capture the asteroid undergoes torque-free motion and 
in general can tumble, and its angular velocity and mass 
properties are known. The spacecraft is an axisymmetrie 
rigid body. The capture mechanism is also axisymmetrie; the 
axes of symmetry of the spacecraft and capture mechanism 
are collinear. Due to symmetry, effective use of the capture 
mechanism does not require any particular relative orienta- 
tion, about the axis of symmetry, between the spacecraft and 
the asteroid. The spacecraft approaches the asteroid along a 
straight line (line of descent) that is fixed in the asteroid, as 
illustrated in Fig. 1. 
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Fig. 1. Illustration of spacecraft approaching asteroid. 


Here, we take the approach of removing completely the 
relative angular velocity between the spacecraft and the 
asteroid, based on the assumption that this offers the best 
chance to avoid damage or destruction of the mechanism 
when capture takes place. Thus, the objectives of spin-up 
control are to make the axis of symmetry of the spacecraft 
parallel to the line of descent, and to eliminate rotational 
motion of the spacecraft relative to the asteroid about the 
line of descent. 

We assume that, after capture, the asteroid is rigidly 
connected to the spacecraft so that together they form one 
rigid body, which is considered unsymmetric to maintain 
generality. The body must undergo de-spin to eliminate an- 
gular velocity relative to the inertial frame. This is necessary, 
for example, to orient solar arrays properly or to aim thrust 
in the direction needed to return to the Earth-Moon system. 
A controller is designed to be asymptotically stable, and stay 
within specified thruster limits. 

The paper is organized as follows. The controllers that 
align the axis of symmetry of the spacecraft with the line of 
descent, and eliminate relative angular velocity between the 
two bodies, are presented first. Next, an analytical expression 
is presented for the control torque required to maintain 
alignment of the axes of symmetry of the spacecraft and 
the special case of an axisymmetrie asteroid. Finally, the 
despin controller is presented, followed by conclusions and 
references. 

II. Attitude Synchronization 

In this section we present controller designs that facilitate 
the synchronization of the spacecraft axis of symmetry and 
the line of descent. To this end, consider a unit vector a 
fixed in the asteroid. A, parallel to the line of descent, and 
having a direction consistent with the descent. Let bi, t> 2 . 
and b 3 be a set of right-handed, mutually orthogonal unit 
vectors fixed in the spacecraft, B, with b 3 parallel to the 
axis of symmetry, and marking the direction in which the 
spacecraft is to advance toward the asteroid. It is desired 
to design a controller to align b 3 with a, and to eliminate 



all relative rotational motion before capture takes place, 
without requiring any particular relative orientation between 
the spacecraft and the asteroid about b 3 . 


A. Kinematics and Dynamics 


The direction of a with respect to the spacecraft reference 
frame is described with two parameters, w\ and w 2 , obtained 
through stereographic projection as discussed in Ref. [4] and 
illustrated in Fig. 2. If the unit vector a is written as 


a = abi + bh 2 + cb 3 


then wi and w 2 are defined as 


w i = 


b 

1 + c 


—a 

w 2 = — — 
1 + c 


(1) 

( 2 ) 



Fig. 2. Illustration of stereographic projection and parameters w\ and W 2 . 


As illustrated in Fig. 2, when w\ = w 2 = 0, b 3 has the 
same direction as a, implying that the axis of symmetry is 
aligned with the line of descent. Thus, the objective of the 
controller design is to drive both W\ and w 2 to zero, and to 
eliminate all relative angular velocity between B and A. 

The problem of pointing the axis of symmetry, b 3 , of 
an axisymmetric spacecraft in the inertial frame is studied 
in Ref. [5], where kinematic equations are given that relate 
the angular velocity of the spacecraft to the rates of w 3 and 
w 2 , as well as several controllers that either exponentially 
or asymptotically stabilize b 3 . The problem in this paper is 
different in that we are dealing with the orientation of b 3 in 
A, which is a rotating frame. 

Let A u> B denote the angular velocity of the spacecraft 
relative to the asteroid. Since a is fixed in frame A, 


A da B da 
di d t 


A w B 


x a = 0 


(3) 


A' d r ) 

where ,y denotes the derivative of a vector with respect 
to time taken in a reference frame X. In addition, Eq. (1) 
leads to 


B da 

df 


= abi + 6b 2 + cb 3 


It can be shown from Eqs. (2), (3), and (4) that w\ and w 2 
are governed by the following kinematic equations. 

£ 

•w i = w 2 S 3 + S 2 wiw 2 + —(1 + w\ — w\) (5a) 

w 2 = ~wiS 3 + 6 iw\w 2 + —(1 + u>f — wf) (5b) 


where 

Si= A u B -bi (i = 1,2,3) (6) 

are the measure numbers (Ref. [6], p. 3) of a uj b for b 1; 
b 2 , and b 3 . 


The dynamics of the rotational motion of the spacecraft are 
described by Euler’s equations. Let Ii (i = 1,2,3) denote 
the principle moments of inertia of the spacecraft. For an 
axisymmetric body, I\ = I 2 . Let 


uti= N ut B - hi (* = 1,2,3) (7) 


denote the b|-b 2 -b :! measure numbers of the angular veloc- 
ity of spacecraft B relative to an inertial frame N. Then, 
Euler’s equations are written as 


Cj i = aiw 2 w 3 + ui (8a) 

w 2 = — aiut 3 uji + u 2 (8b) 

w 3 = u 3 (8c) 

where = (I 2 — I 3 )/Ii and Ui = Af,//, (i = 1, 2, 3). Here, 
Mi, M 2 , and M 3 are the bi-b 2 -b 3 measure numbers of the 
torque acting on the spacecraft. 

The dynamic equations that govern the relative angular 
velocity measure numbers (i = 1,2,3) can then be 
derived. By definition in Eq. (6), 

Si = ( N oj b - N u> A ) • bj (* = 1,2,3) (9) 

Let 

vJi = N ut A • hi (* = 1,2,3) 

denote the b| -b 2 -b 3 measure numbers of the angular veloc- 
ity of the asteroid relative to the inertial frame. Then, 


Si =Ui- Wi (* = 1 , 2 , 3 ) ( 10 ) 

Eqs. (8) and (10) together yield 

Si=aiS 2 S 3 +v i (11a) 

S 2 = —aiS 3 Si + v 2 (lib) 

63 = v 3 (lie) 


where Vi ( * = 1,2,3) are new control variables that are 
related to the original control variables a, (1 - 1,2,3) by 


vi=ui+ ai(w 2 w 3 + S 2 w 3 + w 2 S 3 ) — w\ (12a) 

v 2 = u 2 — a i(w 3 w 3 + S\W 3 + w\S 3 ) — w 2 (12b) 

v 3 =u 3 - w 3 (12c) 

where Wi (i = 1,2,3) are given by 
N 1 

• _ a / n, ,a r t 

at 

= N a A .b l+ V.(Vxb t ) 

= ("a A + ,v u ,1 x ,v u ,! ).b j (* = 1,2,3) 

(13) 


(4) 



where N a A is the angular acceleration of the asteroid 
relative to the inertial frame. 

The asteroid is assumed to have known ratios of moments 
of inertia (Ref. [7]), and undergo a torque-free motion. Thus, 
with an attitude and angular velocity measurement at any 
time, the attitude, angular velocity, and angular acceleration 
of the asteroid can be determined at any other time from 
dynamical and kinematical equations of motion (Ref. [8], 
Secs. 3.1 and 3.4). 


B. Controller Design 

A close examination of Eqs. (5) and (11) reveals that the 
dynamics of <5 3 is self-contained in that it is not affected by 
other states in the system. Thus, the problem can be separated 
into two subsystems, and controllers can be designed to 
stabilize the subsystems individually. 

Subsystem 1: 


5 = — ( 0,1636 + v (14a) 

c- 6 5 , 

w = —iSsw + - + -w (14b) 


Subsystem 2: 

£3 = v 3 (15) 

where Subsystem 1 consists of Eqs. (5), (11a), and (lib), 
while Subsystem 2 consists of Eq. (11c). Here the notation 
of using complex numbers is adopted from Ref. [5]. Namely, 
i = y/—l, 6 = 61 + iS 2 , w = wi + iw 2 , v = v\ + iv 2 , (•) 
denotes complex conjugate, and | • | denotes the magnitude 
of a complex number. 

In Ref. [5], four controllers with exponential or asymptotic 
behavior are developed to align b 3 with n 3 , a unit vector 
fixed in an inertial frame. No control torque is applied about 
b 3 ; therefore, w 3 is constant during the time control is 
applied in the other two directions. Subsystem 1 resembles 
the system studied in Ref. [5]; in general, however, 63 
(the counterpart to constant w 3 ) is not constant. The four 
controllers developed in Ref. [5] are presented in Eqs. 
( 1 6)— ( 1 9) with UJ 3 replaced by 63 . It can be shown in a 
straightforward way that proofs of stability given in Ref. 
[5] hold for <5 3 / 0, and therefore remain applicable to 
Eqs. (16)— (19). For completeness, the four controllers for 
Subsystem 1 are presented as follows. 

Controller 1: 

s s 

V = i(n 6 3 6 + ki{i 6 3 6 - - - -w 2 ) ^ 

— k 2 {5 + k\w) 

with fci > 0 and k 2 > 0 , globally asymptotically stabilizes 
Subsystem 1. 

Controller 2: 


v = iai 6 3 6 - ki (-iS 3 w + ^ + ^ w 2 ) 

- k 2 (S + kiw) - w{l + |w| 2 ) 


with k\ > 0 and k 2 > 0, globally exponentially stabi- 
lizes Subsystem 1 with rate of decay j3/2, where /3 = 
min{fci, 2 ^ 2 }- 


Controller 3: 

v = — kiS — k 2 w (18) 

with ki > 0 and k 2 > 0 , globally asymptotically stabilizes 
Subsystem 1. 

Controller 4: 

v = —k \6 — k 2 {\ + \w\ 2 )w (19) 

with k\ > 0 and k 2 > 0 , globally asymptotically stabilizes 
Subsystem 1. 

The original control variables u± and u 2 can be obtained 
through Eqs. (12a) and (12b) as 

u 1 = i>i — ai(w 2 tU 3 + 6 2 W 3 + W 263 ) + w\ ( 20 a) 

u 2 = v 2 + a\{w\VD 3 + 5 \W 3 + vj\ 63 ) + w 2 ( 20 b) 

Exponential or asymptotic stabilization of Subsystem 1 
ensures that b 3 becomes aligned with a as time progresses. 
These controllers can be readily used in applications such 
as giving b 3 the same direction as a specified vector fixed 
in another rotating satellite, while allowing relative angular 
speed 4 :i to vary freely. However, for retrieving an asteroid, 
it is desirable to drive 63 to zero to avoid damaging the 
capture mechanism. Because of the simple form of Eq. (15) 
associated with Subsystem 2, it is clear that 

V 3 = -k 3 S3 ( 21 ) 

with /,: 3 > 0, globally exponentially stabilizes Subsystem 2 
with rate of decay fc 3 . From Eq. (12c), the original control 
variable u 3 can be obtained as 

U 3 = -k 3 5 3 + TU 3 (22) 

Remark: The controllers presented in this section are 
effective regardless of the initial angle between b 3 and a. 
However, in practice, in the cases where reorientation of b 3 
can be performed without expending much propellant, it can 
be preferable to reorient b 3 first, spin up the spacecraft about 
b 3 , and then apply the feedback controllers. This is because 
the controllers are not designed to minimize propellant usage. 
For example, when the spacecraft is initially at rest, a slow 
reorientation of b 3 can be achieved through brief thruster fir- 
ings, with minimal propellant consumption, at the beginning 
and end of the maneuver. Additionally, the spacecraft can be 
spun up about b 3 once the reorientation is completed. This 
way, the simple attitude maneuvers are separated from more 
complex ones that require feedback control. 

III. CO-SPINNING AXISYMMETRIC SPACECRAFT AND 

Axisymmetric Asteroid 

As the axis of symmetry of the spacecraft asymptotically 
approaches the line of descent, torque must be applied 
continuously to keep the two lines collinear. An analytical 
expression for the torque is available for the special case of 
an axisymmetric asteroid and a line of descent parallel to 
the asteroid’s axis of symmetry. During torque-free motion, 
the central angular momentum H of the asteroid A in a 
Newtonian reference frame N has constant magnitude and 
constant direction in N . The angle </> between H and the 
asteroid’s axis of symmetry is constant. A reference frame C 




Fig. 3. Axisymmetric spacecraft, B, duplicating torque-free rotational 
motion of axisymmetric asteroid. 


(see Fig. 3) can be introduced such that H and the asteroid’s 
axis of symmetry are fixed in C. The angular velocity of 
A in N can be described simply in terms of spin about the 
axis of symmetry, together with precession about the angular 
momentum vector (see Ref. [8], Sec. 3.1). With unit vector £3 
chosen such that it is fixed in C and parallel to the asteroid’s 
axis of symmetry, the angular velocity c u> ' of A in C can 
be written as 

C u A =sc 3 (23) 

where spin speed s is a constant. Moreover, the angular 

velocity N u> G of C in N is given by 

N u c =ph (24) 

where precession speed p is a non-negative constant, and h 
is a unit vector having the same direction as H. Hence, the 
angular velocity of A in N can be expressed simply as 

N ijj a = n uj ° + c u; A = ph + sc 3 (25) 

and the angular acceleration N a A of A in N is given by 

N a A = ^(ph + sc 3 ) =0+s N u G x £3 (26) 

= sph x c 3 = sp sin <f> fi 2 

where £ 2 is a unit vector fixed in C and normal to the plane 
containing H and £3 (see Fig. 3). If the axes of symmetry 
of A and B are to be parallel, then £3 = b 3 , and the central 
inertia dyadic of B can be expressed as 

I = / B (bibi+b 2 b2)TJBb 3 b 3 = I b {ciC 1 +C2C2)+J b c 3 c 3 

(27) 

where Jg and / g are central principal moments of inertia 
of B for, respectively, a line parallel to b 3 , and any line 
perpendicular to b 3 . If B is to duplicate the asteroid’s angular 
velocity and angular acceleration, then N c o B = N us ■ ' and 
N a B = N a ' , and the required control torque is given by 

T = I • N cx A + N u> A xl- n uj a 

= psm<j) [(J B - I B )pcos(j> + J B s] c 2 (28) 


The magnitude of T is constant. The relationship 

cos <t> = Js (29) 

( I ~ J)P 

where I and J are central principal moments of inertia of the 
asteroid, can be used to show that T vanishes if the ratio of 
moments of inertia of the spacecraft happens to be the same 
as that of the asteroid; that is, if /g/Jg = I/J. Otherwise, 
T is the control torque that must be applied to maintain a 
relative angular velocity of zero until capture takes place. 

IV. Despin 

Once the asteroid payload has been securely captured by 
the spacecraft, the composite body (spacecraft and payload) 
undergoes de-spin to eliminate angular velocity relative to 
the inertial frame. For convenience B will denote the newly 
formed rigid body. An asymptotically stable controller is 
designed for applying torque to drive N u> B to 0. No 
requirements are placed on the orientation of B in N at the 
time N u> B reaches 0. It is assumed that the required ori- 
entation can subsequently be obtained by briefly firing RCS 
thrusters at the beginning and end of a slow reorientation 
that consumes a negligible amount of propellant. 

Let I-i (1 =1.2. 3) denote the principle moments of inertia 
of B, and let u>i ( i = 1,2,3) be defined as in Eqs. (7). 
Assume that n RCS thrusters are mounted to the spacecraft, 
and each is capable of producing a moment about the center 
of mass of B. Let F = [Fi, F 2 , . . . , F n \ be a 3 x n matrix, 
the columns of which are the maximum moment produced by 
each thruster. Assume the thrusters can be throttled between 
0 and 1, and let 77 = [771 , 772, • • • , r] n ] T denote the nxl matrix 
of the thruster throttles. In addition, it is assumed that the 
RCS thrusters are configured in such a way that when all 
thrusters are firing at full throttle, the net moment equals to 
zero; i.e., 

FI = 0 (30) 

where 1 = [1, 1, ... , 1] T is the nxl array with all elements 
equal to 1. This is a typical RCS thruster configuration in 
practice. For example, RCS thrusters are often arranged in 
pairs, each of which applies a force couple. 

Euler’s dynamical equations of rotational motion for rigid 
body B can be written as 

Ilo = — u x Iu) + Fri (31) 

where I is a diagonal inertia matrix of B, to = [oji, 0 J 2 , w 3 ] T , 
and L>j x is a skew symmetric matrix used to implement a 
cross product of two vectors. 

Let a Lyapunov function be defined as the rotational 
kinetic energy K of B in N 

K = 1 w T /w (32) 

2 

Hence, the time rate of change of I\ is given by 

K = uj t Ilo = w T (— lo x Iuj + Frj) = uj t Fi 7 (33) 

We seek a strategy for applying control torque that will 
completely dissipate rotational kinetic energy; that is, we 
wish to make K < 0, and K = 0 if and only if uj = 0. 



A fast strategy would be a bang-bang profile for throttle. 
That is. 


Vi = 


if co T Fi < 0 
if up- Fi > 0 


(* = l,2,...,n) (34) 


This strategy, though effective in rapidly dissipating K, 
incurs tremendous chattering especially when uj approaches 
0. Thus, a smooth controller is desired. To this end, we seek a 
controller that not only dissipates K, but also applies torque 
sparingly so as to avoid excessive expenditure of propellant. 
These considerations lead us to an objective function J 
defined as 

J = uj t Ft] + -t] T Rr] (35) 

which will be minimized at each instant of time. Here R > 0 
is an n x n positive definite matrix. The optimal matrix ?/ can 
be found by setting the partial derivative of J with respect 
to ?/ equal to zero; i.e., 

d J r r t 

— = up F + r, T R = 0 (36) 

orj 

which leads to 

r) = —R~ 1 F t uj (37) 


Substituting 77 into Eq. (33) yields 

K = -uj t FR~ 1 F t uj < 0 (38) 


and k = 0 if and only if uj 0. Thus, according 
to Lyapunov theory (Ref. [9]) the controller in Eq. (37) 
asymptotically stabilizes the angular velocity, and minimizes 
the control energy simultaneously. 

Recall that the thruster throttles should be bounded be- 
tween 0 and 1; however, this is not guaranteed in Eq. (37). 
We now proceed to modify Eq. (37) to satisfy the throttle 
constraint and maintain asymptotic stability. To this end, let 
??min be the minimum of all throttles calculated by Eq. (37). 
If 7/m in < 0 , then we modify rj as 

V = -R~ 1 F t uj - 7/ min 1 (39) 


Now, let 7/ max be the maximum of all throttles as given in 
Eq. (39), and if 77 max > 1, then we modify 77 by scaling it 
by 1 / 7 /max ; i-e., 

-R~ 1 F T UJ - 7/min 1 

r) = (40) 

hmax 


It can be deduced that the throttles given in Eq. (40) are all 
between 0 and 1 by construction. Substitute this controller 
into Eq. (33), and we get 


r > T v ( ~ R lFTu ~ ^minl) 
/v — uj r 

?7ma x 

= — uj t FR~ 1 F t uj 

?7max 

< 0 


(41) 


and K = 0 if and only if uo = 0. Here, the assumption in Eq. 
(30) is used. Thus, the controller in Eq. (40) asymptotically 
stabilizes the angular velocity. 

Remark 1: As 7 / m i n is subtracted from each throttle in Eq. 
(39), some thrusters may be fired unnecessarily, resulting 


TABLE I 

RCS Thruster Configuration 


Pod 

Position (m) 

Thrust Directions 

Pi 

r B ' p i = 2bi-2.65b 3 

±b> 2 ? ±£3 

P 2 

r a*p 2 = _ 2 b 1 - 2.65b 3 

±t> 2 , ±63 

P 3 

r s * p 3 = 2 b >2 — 2.65b 3 

±bi, ±63 

Pa 

r s * p i = — 2 t >2 - 2.65b 3 

d=bi, ±t >3 


in needless propellant consumption. The thruster firings 
commanded by Eq. (39) that produce no resultant torques 
need to be eliminated. For example, if thrusters i and j 
with the same maximum thrust produce moments in opposite 
directions, and they are commanded to throttle Vi > Vj > 0, 
then both thrusters are fired unnecessarily. These throttles can 
be clipped by setting the new r]i to r/> — rjj, and the new r/j 
to 0. Doing so does not change the commanded torque, and 
thus does not affect the stability of the closed-loop system. 

Remark 2: The special case where the thrusters are 
arranged in collocated pairs, each of which provides thrust in 
opposite directions, can be equivalently described as having 
n/2 thrusters with throttles between —1 and 1. Let rfmax 
denote the maximum absolute value of the throttles given in 
Eq. (37). If > 1, then the controller in Eq. (37) can be 
modified to satisfy the throttle bounds as 


V = 


R~ 1 F t uj 

Vmax 


(42) 


It is straightforward to see that the asymptotic stability of 
the controller is maintained. 


V. Numerical Examples 

In this section, numerical examples will be presented to 
show the effectiveness of the controllers. The spacecraft 
design is adopted from Ref. [1]. The mass of the spacecraft 
is 18,000 kg. The cylindrical bus is axisymmetric, with a 
diameter of 2.7 m and a height of 5.9 m. Four pods, each 
containing four RCS thrusters, are mounted to the bus on 
struts as shown in Fig. 5 of Ref. [1] (p. 23). Each pod has a 
radial distance of 2 m from the cylinder’s axis of symmetry, 
and is 0.3 m from the cylinder’s end. Each thruster can apply 
200 N of force, and has a specific impulse Isp = 287 s. The 
position vector r B Pi from the mass center of the bus, B* , 
to each of the four pods Pi ( i - 1 , 2 , 3, 4), and directions of 
thrust, are reported in Table 1. 

A. Example 1: Attitude Synchronization 

Use of the feedback controller to align b 3 with a is 
demonstrated in an example involving a tumbling unsym- 
metric asteroid, taken to be an ellipsoid with semidiameters 
a = 5 m, /? = 3 m, and 7 = 4 m, and a uniform 
density of 2000 kg/m 3 . The angular velocity of A in N 
is characterized at / = 0 by the values N u> A • ax = 0.6 
deg/s, N u> A • £12 = 0.6 deg/s, and N u> A • a 3 = 6.0 deg/s 
(60 revolutions per hour). The line of descent, a, is chosen 
to have the same direction as a 3 . 

Although four controllers are presented in Section II-B, 
only the linear controller in Eq. (18) is demonstrated here. 



It is assumed here that the procedure of reorienting the 
spacecraft and spinning up about the axis of symmetry is 
executed before t = 0. That is, at t = 0, b :{ is aligned 
with a. 3 , and the spacecraft spins about b :i at 6.0 deg/s. In 
addition, the angle between bi and ai is arbitrarily chosen 
to be 60 deg. The controller consisting of Eqs. (18), (20), 
and ( 22 ) is applied at t = 0 . 

Results are shown in Fig. 4. Figure 4(a) shows the measure 
numbers Si ( i = 1,2,3) of the angular velocity of the 
spacecraft relative to the asteroid. The controller is able to 
drive the errors to zero in about 200 s. The angle between the 
spacecraft axis of symmetry and the line of descent, which 
can be expressed as cos _ 1 (b 3 • £ 13 ), is shown in Fig. 4(b). 
The angle is driven to zero within about 200 s. The torque 
commanded by the controller is shown in Fig. 4(c). Notice 
that after 200 s torque is required continuously to keep the 
angular velocity errors at zero, and S 3 aligned with as. 


(a) 



Time (s) 


Fig. 4. (a) Relative angular velocity measure numbers 4, = A u> B • b, 

(i = 1,2,3), (b) angle between the axis of symmetry and the line of descent, 
(c) control torque required for attitude synchronization. 

B. Example 2: Despin 

The despin controller is demonstrated with the spacecraft 
and asteroid introduced in Example 1. The result of capture is 
that b 3 has the same direction as a 3 , the top of the spacecraft 
is in contact with the asteroid, and a 60-deg relative angle 
about b 3 is present, as discussed in Example 1. Unit vectors 
bi, t> 2 , and t >3 fixed in the composite body B are chosen to 
be parallel to the diameters of the ellipsoidal asteroid. The 
angular velocity of B in N is characterized at t - 0 by the 
values N uj B • h, (i = 1 , 2 , 3) equal to 0.6 deg/s, 0.6 deg/s, 
and 6.0 deg/s, respectively. 

Because the mass of the captured asteroid is significantly 
greater than that of the spacecraft, the center of mass of B 
is close to that of the asteroid. In this case, moment in the 
direction of bi or b 2 is applied with an arm that is 4.2 times 
as large as the arm associated with moment in the direction 
of t> 3 . In Ref. [2], tumbling is intentionally induced in order 
to take advantage of larger moment arms. In the present 
example, the largest component of angular velocity is in the 



“0 500 1000 1500 

Time (s) 

Fig. 5. Angular velocity of the composite body during despin. 

direction of b 3 , which is parallel to the intermediate axis 
of inertia, and the initial conditions correspond to a state of 
tumbling. By waiting for 251.7 s until the angular momentum 
is transferred to the plane spanned by bi and t> 2 , and then 
applying the despin controller, one can take advantage of 
the larger moment arms. The controller drives the angular 
velocity asymptotically to zero, as shown in Fig. 5, and 
34 kg of propellant in consumed. The cost of not waiting 
for angular momentum transfer is 64 kg of propellant, an 
increase of 88 %. 

VI. Conclusion 

Attitude control associated with capturing a free-flying 
asteroid using an axisymmetric spacecraft is considered. It 
has been shown that the axis of symmetry of the spacecraft 
can be aligned asymptotically with a line of descent fixed 
in an asteroid, and the relative angular velocity between the 
two bodies can be made to vanish. An analytical expression 
is obtained for the torque required to keep the spacecraft axis 
of symmetry aligned with that of an axisymmetric asteroid. 
After the asteroid is securely captured, the angular velocity of 
the composite body relative to an inertial frame can be driven 
asymptotically to zero. In some cases propellant consumption 
can be reduced by waiting for angular momentum transfer 
to occur before applying the despin controller, because then 
advantage can be taken of longer momentum arms about 
certain axes. 
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